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We present randomized approximation algorithms for multi-criteria Max-TSP. 
For Max-STSP with k > 1 objective functions, we obtain an approximation ratio 
of t — £ for arbitrarily small e > 0. For Max-ATSP with k objective functions, we 
obtain an approximation ratio of — £ - 

1 Multi-Criteria Traveling Salesman Problem 

1.1 Traveling Salesman Problem 

The traveling salesman problem (TSP) is one of the most fundamental problems in combina- 
torial optimization. Given a graph, the goal is to find a Hamiltonian cycle of minimum or 
maximum weight. We consider finding Hamiltonian cycles of maximum weight (Max-TSP). 

An instance of Max-TSP is a complete graph G = (V, E) with edge weights w : E — » N. The 
goal is to find a Hamiltonian cycle of maximum weight. The weight of a Hamiltonian cycle 
(or, more general, of a subset of E) is the sum of the weights of its edges. If G is undirected, 
we speak of Max-STSP (symmetric TSP). If G is directed, we have Max-ATSP (asymmetric 
TSP). 

Both Max-STSP and Max-ATSP are NP-hard and APX-hard. Thus, we are in need of 
approximation algorithms. The currently best approximation algorithms for Max-STSP and 
Max-ATSP achieve approximation ratios of 61/81 and 2/3, respectively [2,5]. 

Cycle covers are an important tool for designing approximation algorithms for the TSP. 
A cycle cover of a graph is a set of vertex-disjoint cycles such that every vertex is part of 
exactly one cycle. Hamiltonian cycles are special cases of cycle covers that consist of just one 
cycle. Thus, the weight of a maximum-weight cycle cover is an upper bound for the weight 
of a maximum-weight Hamiltonian cycle. In contrast to Hamiltonian cycles, cycle covers of 
minimum or maximum weight can be computed efficiently using matching algorithms [1]. 

1.2 Multi-Criteria Optimization 

In many optimization problems, there is more than one objective function. Consider buying 
a car: We might want to buy a cheap, fast car with a good gas mileage. How do we decide 
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which car suits us best? With multiple criteria involved, there is no natural notion of a best 
choice. Instead, we have to be content with a trade-off. The aim of multi-criteria optimization 
is to cope with this problem. To transfer the concept of an optimal solution to multi-criteria 
optimization problems, the notion of Pareto curves was introduced (cf. Ehrgott [3]). A Pareto 
curve is a set of solutions that can be considered optimal. 

More formally, a /c-criteria optimization problem consists of instances /, solutions sol(X) for 
every instance X £ I, and k objective functions w±, . . . ,Wk that map X £ I and Y € sol(X) 
to N. Throughout this paper, our aim is to maximize the objective functions. We say that 
a solution Y E sol(A) dominates another solution Z € sol(A) if Wi(Y, X) > Wi(Z,X) for all 
i 6 [k] = {1, . . . , k} and w i(Y, X) > Wi(Z, X) for at least one i. This means that Y is strictly 
preferable to Z. A Pareto curve (also known as Pareto set or efficient set) for an instance 
contains all solutions of that instance that are not dominated by another solution. 

Unfortunately, Pareto curves cannot be computed efficiently in many cases: First, they 
are often of exponential size. Second, because of straightforward reductions from knapsack 
problems, they are NP-hard to compute even for otherwise easy problems. Thus, we have to 
be content with approximate Pareto curves. 

For simpler notation, let w(Y, X) = (wi(Y, X), . . . ,Wk(Y,X)). We will omit the instance X 
if it is clear from the context. Inequalities are meant component-wise. A set V Q sol(A) of 
solutions is called an a approximate Pareto curve for X € / if the following holds: For every 
solution Z G sol(A), there exists a Y € V with w(Y) > aw(Z). We have a < 1, and a 1 
approximate Pareto curve is a Pareto curve. (This is not precisely true if there are several 
solutions whose objective values agree. However, in our case this is inconsequential, and we 
will not elaborate on this for the sake of clarity.) An algorithm is called an a approximation 
algorithm if, given the instance X, it computes an a approximate Pareto curve. It is called a 
randomized a approximation algorithm if its success probability is at least 1/2. This success 
probability can be amplified to 1 — 2~ m by executing the algorithm m times and taking the 
union of all sets of solutions. (We can also remove solutions from this union that are dominated 
by other solutions in the union, but this is not required by the definition of an approximate 
Pareto curve.) 

Papadimitriou and Yannakakis [10] showed that (1 — e) approximate Pareto curves of size 
polynomial in the instance size and 1/e exist. The technical requirement for the existence is 
that the objective values of solutions in sol(A) are bounded from above by 2 p ( N ^ for some 
polynomial p, where N is the size of X. This is fulfilled in most natural optimization problems 
and in particular in our case. 

A fully polynomial time approximation scheme (FPTAS) for a multi-criteria optimization 
problem computes (1 — e) approximate Pareto curves in time polynomial in the size of the 
instance and 1/e for all e > 0. Papadimitriou and Yannakakis [10], based on a result of 
Mulmuley et al. [9] , showed that multi-criteria minimum- weight matching admits a randomized 
FPTAS, i. e., the algorithm succeeds in computing a (1 — e) approximate Pareto curve with 
constant probability. This randomized FPTAS yields also a randomized FPTAS for the multi- 
criteria maximum- weight cycle cover problem [8], which we will use in the following. 

Manthey and Ram [6, 8] designed randomized approximation algorithms for several variants 
of multi-criteria Min-TSP. However, they leave it as an open problem to design any approxi- 
mation algorithm for Max-TSP. 
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1.3 New Results 



We devise the first approximation algorithm for multi-criteria Max-TSP. For A;-criteria Max- 
STSP, we achieve an approximation ratio of \ — e for arbitrarily small e > 0. For /c-criteria 
Max-ATSP, we achieve — e. Our algorithm is randomized. Its running-time is polynomial 
in the input size and 1/e and exponential in the number k of criteria. However, the number of 
different objective functions is usually a small constant. 

The main ingredient for our algorithm is a decomposition technique for cycle covers and a 
reduction from /c-criteria instances to {k — l)-criteria instances. 

2 Outline and Idea 

A straight-forward 1/2 approximation for mono-criterion Max-ATSP is the following: First, 
we compute a maximum-weight cycle cover C. Then we remove the lightest edge of each cycle, 
thus losing at most half of C's weight. In this way, we obtain a collection of paths. Finally, we 
add edges to connect the paths to get a Hamiltonian cycle. For Max-STSP, the same approach 
yields a 2/3 approximation since the length of every cycle is at least three. 

Unfortunately, this does not generalize to multi-criteria Max-TSP for which "lightest edge" is 
usually not well defined: If we break an edge that has little weight with respect to one objective, 
we might lose a lot of weight with respect to another objective. Based on this observation, the 
basic idea behind our algorithm and its analysis is the following case distinction: 

Light-weight edges: If all edges of our cycle cover contribute only little to its weight, then 
removing one edge does not decrease the overall weight by too much. Now we choose the 
edges to be removed such that no objective loses too much of its weight. 

Heavy-weight edges: If there is one edge that is very heavy with respect to at least one objective, 
then we take only this edge from the cycle cover. In this way, we have enough weight for 
one objective, and we proceed recursively on the remaining graph with k — 1 objectives. 

In this way, the approximation ratio for /c-criteria Max-TSP depends on two questions: 
First, how well can we decompose a cycle cover consisting solely of light-weight edges? Second, 
how well can (k — l)-criteria Max-TSP be approximated? We deal with the first question 
in Section 3. In Section 4, we present and analyze our approximation algorithms, which 
also gives an answer to the second question. Finally, we give evidence that the analysis of the 
approximation ratios is tight and point out some ideas that might lead to better approximation 
ratios (Section 5). 

3 Decompositions 

Let a <G (0,1], and let C be a cycle cover. We call a collection P C C of paths an a- 
decomposition of C if w(P) > aw(C). (Remember that all inequalities are meant component- 
wise.) In the following, our aim is to find a-decompositions of cycle covers consisting solely of 
light-weight edges, that is, w{e) < aw{C) for all e G C. 

Of course, not every cycle cover possesses an a-decomposition for every a. For instance, a 
single directed cycle of length two, where each edge has a weight of 1 shows that a = 1/2 is 
best possible for a single objective function in directed graphs. On the other hand, by removing 
the lightest edge of every cycle, we obtain a 1/2-decomposition. 
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For undirected graphs and k = 1, a = 2/3 is optimal: We can find a 2/3-decomposition by 
removing the lightest edge of every cycle, and a single cycle of length three, where each edge 
weight is 1, shows that this is tight. 

More general, we define af G (0, 1] to be the maximum number such that every directed 
cycle cover C with w(e) < af-w(C) for all e G C possesses an a^-decomposition. Analogously, 
G (0, 1] is the maximum number such that every undirected cycle cover C with w(e) < 
a\ • w{C) possesses an a^-decomposition. We have af = \ and a\ = |, as we have already 
argued above. We also have a\ > af and a\ < ajj^ as well as af < af_ 1 - 

3.1 Existence of Decompositions 

In this section, we investigate for which values of a such a-decompositions exist. In the 
subsequent section, we show how to actually find good decompositions. We have already dealt 
with a" and af . Thus, k > 2 remains to be considered in the following theorems. In particular, 
only k > 2 is needed for the analysis of our algorithms. 

Let us first normalize our cycle covers to make the proofs in the following a bit easier. For 
directed cycle covers C, we can restrict ourselves to cycles of length two: If we have a cycle c 
of length £ with edges ei, . . . , eg, we replace it by [£/2\ cycles faj-i, for j = 1, . . . , \£/2\. 
If I is odd, then we add a edge ee + i with w(ee + i) = and add the cycle (eg, e^+i). (Strictly 
speaking, edges are 2-tuples of vertices, and we cannot simply reconnect them. What we mean 
is that we remove the edges of the cycle and create new edges with the same names and weights 
together with appropriate new vertices.) We do this for all cycles of length at least three and 
call the resulting cycle cover C . Now any a-decomposition P' of the new cycle cover C yields 
an a-decomposition P of the original cycle cover C by removing the newly added edges e^+i: 
In C, we have to remove at least one edge of the cycle c to obtain a decomposition. In C", we 
have to remove at least [£/2\ edges of c, thus at least one. Furthermore, if w{e) < a ■ w(C) for 
every e G C, then also w(e) < a ■ w(C) for every e G C since we kept all edge weights. This 
also shows w(P) = w(P r ). 

We are interested in a-decompositions that work for all cycle covers with k objective func- 
tions. Thus in particular, we have to be able to decompose C. The consequence is that if every 
directed cycle cover that consists solely of cycles of length two possesses an a-decomposition, 
then every directed cycle cover does so. 

For undirected cycle covers, we can restrict ourselves to cycles of length three: We replace a 
cycle c = (e±, . . . , ei) by [£/3\ cycles (63^-2, £3j-i, for 1 < j < [£/3\. If £ is not divisible by 
three, then we add one or two edges e^+i, e^ + 2 to form a cycle of length three with the remaining 
edge(s). Again, every a-decomposition of the new cycle cover yields an a-decomposition of the 
original cycle cover. 

In the remainder of this section, we assume that all directed cycle covers consist solely of 
cycles of length two and all undirected cycle covers consist solely of cycles of length three. 
Both theorems are proved using the probabilistic method. 

3.1.1 Undirected Cycle Covers 

For the proof of Theorem 3.2 below, we use Hoeffding's inequality [4, Theorem 2], which we 
state here in a slightly modified version. 
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Lemma 3.1 (Hoeffding's inequality). Let X\, . . . ,X n be independent random variables, where 
Xj assumes values in [aj,bj]. Let X = Y^=i-^j- Then 

F(X < E(X) -t)< exp 

Theorem 3.2. For all k > 2, we have a% > ^. 

Proof. Let C be any cycle cover and u>i , . . . , Wk be k objective functions. First, we scale the 
edge weight such that Wi(C) = k for all i. Thus, Wi(e) < 1 for all edges e of C since the 
weight of any edge is at most a 1/k fraction of the total weight. Second, we can assume that 
C consists solely of cycles of length three. 

Let ci, . . . , c m be the cycles of C and let ej, e|, e| be the three edges of Cj. We perform the 
following random experiment: We remove one edge of every cycle independently and uniformly 
at random to obtain a decomposition P. Fix any i € [k]. Let Xj be the weight with respect 
to Wi of the path in P that consists of the two edges of Cj. Then E(Xj) = 2wi(cj)/3. Let 
X = T.J=i x j- Th en E( Wi (X)) = 2 Wi (C)/3 = 2fc/3. 

Every Xj assumes values between aj = min{u>j (ej ) +iWj(ej), itfj(ej) + u?j(e|), Wj(e|) + u>j(e|)} 
and bj = max{tt>i(ej) + u>i(e|), Wj(ej) + Wj(e|), Wj(e|) + Wi(ej)}. Since the weight of each edge 
is at most 1, we have bj — aj < 1. Since the sum of all edge weights is A:, we have 

mm m 

k>^bj>^bj- aj > J2( b i ~ a if- 
j=i j=i j=i 

Let us estimate the probability of the event that X < 1, which corresponds to Wi(P) < 1. 
If F(X < 1) < 1/k, then, by a union bound, we have F(3i : Wi(P) < 1) < 1. Thus, P(Vi : 
Wi{P) > 1) > 0, which implies the existence of a l//c-decomposition. By Hoeffding's inequality, 

( 2fc (2k \\ ( 2(f-l) 2 \ 

F(X<l)=F(x< T -( y -l))<exp^AA_Xj =:pk . 

We have ^4 « 0.2494, ^5 ss 0.11, and p$ w 0.05. Thus, for = 4,5,6, and also for all larger 
values of k, we have p& < 1/k, which implies the existence of a l/fc-decomposition for k > 4. 
The cases k = 2 and /c = 3 remain to be considered since ~ 0.51 > 1/3 and p2 ~ 0.89 > 1/2. 
The bound for a^, follows from Lemma 3.3 below, which does not require Wi(e) < • Wi(C). 

Let us show > 1/3. This is done in a constructive way. First, we choose from every 
cycle Cj the edge that maximizes w% and put it into P'. The set P' will become a subset 
of P. Then w^{P') > 1. But we can also have some weight with respect to w\ or W2- Let 
Si = wi(P') and S2 = W2(P')- If 6i > 1, then wi does not need any further attention. 

Let C = C\ P' . We have Wi(C) = 3 — Si for i = 1, 2, and C" consists solely of paths of 
length two. Of every such path, we can choose at most one edge for inclusion in P. (Choosing 
both would create a cycle.) Let ej,e| be the two edges of Cj with W2(eJ) > ^(ej)- Now we 
proceed by considering only W2- Let Q,Q' be initially empty sets. For all j = 1, ... ,m, if 
wi{Q) > W2{Q'), then we put (the heavier edge) e 2 j into Q' and e] into Q. If W2{Q) < W2{Q'), 
then we put e| into Q and ej into Q'. 

Both P' U Q and P' U Q' are decompositions of C. We claim that at least one has a 
weight of at least 1 with respect to all three objectives. Since ws(P') > 1, this holds for both 
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with respect to W3. Furthermore, \w2(Q) — w 2 (Q')\ < 1 since W2(e) < 1 for all edges. We 
have w 2 (Q) + w 2 (Q') = 3 - 5 2 . Thus, mm{w 2 (Q),w 2 (Q')} > ^ - \ > 1 - f ■ This implies 
w 2 (P'uQ) > 1 and w 2 (P'uQ') > 1. Hence, with respect to w 2 and w 3 , both P'UQ and P'UQ 
will do. The first objective w\ remains to be considered. We have max{wi(Q), wi(Q')} > 
Choosing either P = P' U Q or P = P' U Q' results in wi(P) > 81 + ^ > 1. □ 

For undirected graphs and k = 2, we do not need the assumption that the weight of each 
edge is at most times the weight of the cycle cover. Lemma 3.3 below immediately yields a 
(1/2 — e) approximation for bi-criteria Max-STSP: First, we compute a Pareto curve of cycle 
covers. Second, we decompose each cycle cover to obtain a collection of paths, which we then 
connect to form Hamiltonian cycles. The following lemma can also be generalized to arbitrary 
k (Lemma 3.6). 

Lemma 3.3. For every undirected cycle cover C with edge weights w = (wi,w 2 ), there exists 
a collection P C C of paths with w(P) > w(C)/2. 

Proof. Let c be a cycle of C consisting of edges e±, e 2 , e%. Since we have three edges, there exists 
one edge ej that is neither the maximum-weight edge with respect to w\ nor the maximum- 
weight edge with respect to w 2 . We remove this edge. Thus, we have removed at most half 
of c's weight with respect to either objective. Consequently, we have kept at least half of c's 
weight, which proves a 2 > 1/2. □ 

3.1.2 Directed Cycle Covers 

For directed cycle covers, our aim is again to show that the probability of having not enough 
weight in one component is less than 1/k. Hoeffding's inequality works only for k > 7. We use 
a different approach, which immediately gives us the desired result for k > 6, and which can 
be tweaked to work also for small k. 

Theorem 3.4. For all k > 2, we have af > jrpy. 

Proof. As argued above, we can restrict ourselves to cycle covers consisting solely of cycles of 
length two. We scale the edge weights to achieve Wi(C) = k + 1 for all i G [A;]. This implies 
wi(e) < 1 for all edges e G C. 

Of every cycle, we randomly choose one of the two edges and put it into P. Fix any i 6 [A;]. 
Our aim is to show that ¥(wi(P) < 1) < 1/A:, which would prove the existence of an at- 
decomposition. Let ci,...,c m be the cycles of C with Cj = {ej,fj)- Let Wi{ej) = aj and 
Wi(fj) = bj. We assume a,j < bj for all j E [m]. Let S = Xlj=i a j- Then, no matter which 
edges we choose, we obtain a weight of at least 5. Hence, if S > 1, we are done. Otherwise, we 
have 5 < 1 and replace bj by bj — aj and aj by 0. Then we only need additional weight 1 — 5, 
and our new goal is to prove F(wi(P) < 1 — 5) < 1/k. 

This boils down to the following random experiment: We have numbers b\, . . . , b m G [0, 1] 
with X^j=i bj = k + 1 — 25. Then we choose a set / C [m] uniformly at random. For such an 
/, we define (by abusing notation) w(I) = J2jei ^ e nave *° snow ¥{w{I) < 1 — 5) < 1/k. 

To this aim, let Ci, . . . ,C Z C [m] with z = be pairwise disjoint sets with w(Cg) G 

[1 — 5,2 — 5). Such sets exist: We select arbitrary elements for C\ until w(C\) G [1 — 5,2 — 5). 
This can always be done since bj < 1 for all j. Then we continue with C 2 , C3, and so on. If 
we have already z — 1 such sets, then 

w(d U . . . U C*_i) < (2 - 5) ■ (z - 1) < (2 - 8) ■ I < jfe - 6 
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since k > 2. Thus, at least weight k + 1 — 25 — (k — 5) = 1 — 5 is left, which suffices for C z . 

The sets C±, . . . , C z do not necessarily form a partition of [m]. Let C" = [to,] \ [C\ U . . . U C z ). 
We will have to consider C once in the end of the proof. 

Now consider any I,JC [to] . We say that / ~ J if 

I = JAC ei AC e2 A...AC ey 

for some C^, . . . ,Ci y . Here, A denotes the symmetric difference of sets. The relation ~ is 
an equivalence relation that partitions all subsets of [to] into 2 m ~ z equivalence classes, each of 
cardinality 2 Z . Let [I] = {J C [to] | J ~ /}. 

Lemma 3.5. For every I C [to], there are at most two sets J G [I] with w(J) < 1 — 5. 

Proof. Without loss of generality assume that w(I) = minj e ^] w(J). If w(I) > 1 — 5, then 
there is nothing to show. Otherwise, consider any J = IAC^A . . . ACe y € [I] with y > 2: 

y y y 

w(J)>^Tw(c £p \i)> ^Tw(c £p ) -^ w ( Cep ni)>i-5. 
p=i p=i p=i 

V v ' V v ' 

>y-(l— <5)>2— 2<5 <w(I)<l-S 

We conclude that the only possibility for other sets J € [I] with w(J) < 1 — 5 is J = IACn for 
some I. We prove that there is at most one such set by contradiction. So assume that there 
are J x = I Ad and J 2 = IAC 2 with w{Ji),w{J 2 ) < 1 - 5. Then w(Ji) > w{d) - w(d n 
J) + w(C 2 n /) and w(J 2 ) > w{C 2 ) - w{d n I) + w{d n I). Thus, 

2-25> w(Ji) + w{J 2 ) > w{d) + w(C 2 ) >2-25, 

a contradiction. □ 

A consequence of Lemma 3.5 is F(w(I) < 1 — 5) < 2~ z+1 = 2 T ^— 1 _ This is less than 1/k 
for k > 6. The cases k G {2, 3, 4, 5} need special treatment. 

Let us treat k G {2,4} first. Here 2 T ^ 1 +i = jy^ which is almost good enough. To prove 
F(w(I) < 1 — 5) < 1/k, we only have to find a set / such that at most one set J G [I] has 
w(J) < 1 — 5. We claim that is such a set: Of course w($) = < 1 — 5. But for any other 
J G [0], we have 

J = 0AQ 1 A...AQ y =C £l U...Ud y 

for some Q 15 . . . , C^. The latter equality holds since the sets C±, . . . , C z are disjoint. Thus 
w(J) >y(l-5)>l-5. 

To finish the proof, we consider the case k G {3,5}. For this purpose, we consider / and 
I = [m]\I simultaneously. The classes [I] and [I] are disjoint: [I] = [I] would imply C' = 0. 
Then YL\=\ w (Ce) = k + 1 — 25. Since k is odd, we have z = ^p. Thus, since k + 1 > 4, there 
must exist an £ with 

v ^ - k±i k + l 

which contradicts w (Ce) < 2 — 5. 
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We show that the number of sets J G [I] U [I] with u>(J) < 1 — 5 is at most two. This 
would prove the result for k G {3, 5} since this would improve the bound to ¥(w(I) < 1 — 5) < 

2-(*+i)+i = < l/fc. 

If we had more than two sets J G [I] U [/] with io(J) < 1 — 5, we can assume that we have 
two such sets in [I]. (We cannot have more than two such J due to Lemma 3.5.) We assume 
that these two sets are / and V = IAC\. Now consider any J G [/]. Since k is odd and k + 1 
is even, we have 

z 

w(J) < J^w(C e ) + max{w{I),w(l')} < (2 - 5) ■ 

^L, — ' <i-s 

<{2-S)-(z-l) 

- 1 = k - — ■ 5 < k - 5. 
2 

Thus, all sets J G [I] have a weight of less than k — 5. This implies w(J) = k + 1 — 25 — w(J) > 
1 — 5 for all J € [/]. Thus, if [I] contains two sets whose weight is less than 1 — 5, then [I] 
contains no such set. □ 



k + 1 
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3.1.3 Improvements and Generalizations 

To conclude this section, let us discuss some improvements of the results of this section. First, 
as a generalization of Lemma 3.3, cycle covers without cycles of length at most k can be 1/2- 
decomposed. This, however, does not immediately yield an approximation algorithm since 
finding maximum-weight cycle covers where each cycle must have a length of at least k is NP- 
and APX-hard for k > 3 in directed graphs and for k > 5 [7]. 

Lemma 3.6. Let k > 1, and let C be an arbitrary cycle cover such that the length of every 
cycle is at least k+1. Then there exists a collection PCC of paths with w(P) > w(C)/2. 

Proof. The proof is similar to the proof of Lemma 3.3. Let c be any cycle of C. For each 
j 6 [k], we choose one edge of c that maximizes Wi for inclusion in P. Since c has at least k + 1 
edges, this leaves us (at least) one edge for removal. □ 

Figure 1 shows that Theorems 3.2 and 3.4, respectively, are tight for k = 2. Due to these 
limitations for k = 2, proving larger values for a\ or af does not immediately yield better 
approximation ratios (see Section 5). However, for larger values of k, Hoeffding's inequality 
yields the existence of fi(l/ log /^-decompositions. Together with a different technique for 
heavy-weight cycle covers, this might lead to improved approximation algorithms for larger 
values of k. 

Lemma 3.7. We have a%, af G 0,(1/ log k). 

Proof. Let A = cln k + d for some sufficiently large constants c and d. Since a\ > af, we 
can restrict ourselves to directed graphs. Using the notation of Theorem 3.2, we have to show 
that ¥(X < 1) < l/fc, where X = YOjLi-Xj and Xj assumes values in the interval [aj,bj], 
bj < a,j + 1, Yl^Liibj + a j) 2 < A an d E(-X") = A/2. We use Hoeffding's inequality and plug in 
t = A/2-1: 

P(^<i)<«p(-^)=-p(-4 + 2 -|)<i. 
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(a) a% < 1/3. (b) o% < 1/2- 



Figure 1: Examples that limit the possibility of decomposition. 

P <— Decompose(C, w, k, a) 
input: cycle cover C, edge weights w, k > 2, w(e) < a ■ w(C) for all e G C 
output: a collection P of paths 

1: obtain w' from to by scaling each component such that w'^C) = 1/a for all i 

2: normalize C to C as described in the text such that C consists solely of cycles of length 
three (undirected) or two (directed) 

3: while there are cycles c and d in C' with w'(c) < 1/2 and w'(d) < 1/2 do 

4: combine c and c' to c with w/(c) = u>'(c) + w'(d) 

5: replace c and c' by c in C' 

6: try all possible combinations of decompositions 

7: choose one P' that maximizes min iG [fc] w[(P) 

8: translate P' C C' back to obtain a decomposition P C C 

9: return P 

Algorithm 1: A deterministic algorithm for finding a decomposition. 

□ 

3.2 Finding Decompositions 

While we know that decompositions exist due to the previous section, we have to find them 
efficiently in order to use them in our approximation algorithm. We present a deterministic 
algorithm and a faster randomized algorithm for finding decompositions. 

3.2.1 Deterministic Algorithm 

Decompose (Algorithm 1) is a deterministic algorithm for finding a decomposition. The 
idea behind this algorithm is as follows: First, we scale the weights such that w(C) = 1/a. 
Then w(e) < 1 for all edges e £ C. Second, we normalize all cycle covers such that they 
consist solely of cycles of length two (in case of directed graphs) or three (in case of undirected 
graphs). Third, we combine very light cycles as long as possible. More precisely, if there are 
two cycles c and d such that w'(c) < 1/2 and w'{d) < 1/2, we combine them to one cycle c with 
w'(c) < 1. The requirements for an a-decomposition to exist are still fulfilled. Furthermore, 
any a-decomposition of C immediately yields an a-decomposition of C. 

The proof of the following lemma follows immediately from the existence of decompositions 
(Theorems 3.2 and 3.4). 

Lemma 3.8. Let k > 2. Let C be an undirected cycle cover and w\, . . . ,Wk be edge weights 
such that w(e) < a\ ■ w(C). Then Decompose(C, w, k, a%) returns a collection P of paths 
withw{P)>al-w{C). 
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P <— RandDecompose(C, w, fc, a) 




input: cycle cover C, edge weights w = (w±, . . . , Wk), k > 2, w 


e) < a ■ w(C) for all e e C 


output: a collection P of paths with w(P) > a ■ w{C) 




1 


if A; > 6 then 




2 


repeat 




3 


randomly choose one edge of every cycle of C 




4 


remove the chosen edges to obtain P 




5 


until w(P) > a ■ w{C) 




6 


else 




7 


P <— Decompose(C, w, k, a) 






Algorithm 2: A randomized algorithm for finding 


a decomposition. 



Xei C be a directed cycle cover and w i, . . . , Wk be edge weights such that w(e) < af ■ w(C). 
Then Decompose(C, w, k, af ) returns a collection P of paths with w(P) > af ■ w(C). 

Let us also estimate the running-time of Decompose. The normalization in lines 1 to 5 can 
be implemented to run in linear time. Due to the normalization, the weight of every cycle is 
at least 1/2 with respect to at least one w\. Thus, we have at most 2k /at% cycles in C' in the 
undirected case and at most 2k /oti cycles in C in the directed case. In either case, we have 
0(k 2 ) cycles. All of these cycles are of length two or of length three. Thus, we find an optimal 
decomposition, which in particular is an a\ or o^-decomposition, in time linear in the input 
size and exponential in k. 

3.2.2 Randomized Algorithm 

By exploiting the probabilistic argument of the previous section, we can find a decomposition 
much faster with a randomized algorithm. RandDecompose (Algorithm 2) does this: We 
choose the edges to be deleted uniformly at random for every cycle. The probability that 
we obtain a decomposition as required is positive and bounded from below by a constant. 
Furthermore, as the proofs of Theorems 3.2 and 3.4 show, this probability tends to one as k 
increases. For k > 6, it is at least approximately 0.7 for undirected cycle covers and at least 1/4 
for directed cycle covers. For k < 6, we just use our deterministic algorithm, which has linear 
running-time for constant k. The following lemma follows from the considerations above. 

Lemma 3.9. Let k > 2. Let C be an undirected cycle cover and wi, . . . ,Wk be edge weights 
such that w(e) < a\ ■ w(C). Then RandDecompose(C, w, k, a%) returns a collection P of 
paths with w(P) > a\ ■ w(C). 

Let C be a directed cycle cover and w±, . . . ,Wk be edge weights such that w(e) < af ■ w(C). 
Then RandDecompose(C, w, k, af) returns a collection P of paths with w(P) > af ■ w(C). 

The expected running-time of RandDecompose is 0(\C\). 

4 Approximation Algorithms 

Based on the idea sketched in Section 2, we can now present our approximation algorithms 
for multi-criteria Max-ATSP and Max-STSP. However, in particular for Max-STSP, some 
additional work has to be done if heavy edges are present. 
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Vtsp *- MC-MAxATSP{G,w,k,e) 
input: directed complete graph G = (V,E), k > 1, edge weig hts w : E -> N k , e > 
output: approximate Pareto curve "Ptsp f° r fe-criteria Max-TSP 



1: if k = 1 then 

2: compute a 2/3 approximation "Ptsp 
3: else 

4: compute a (1 — e) approximate Pareto curve C of cycle covers 

5: Vtsp <- 

6: for all cycle covers C G C do 

7: if iw(e) < • w(C) for all edges e <G C then 

8: P <— DECOMPOSE(C, 10, fc) 

9: add edges to P to form a Hamiltonian cycle H 

10: add PT to Ptsp 

11: else 

12: let e = (u, v) G C be an edge with w(e) j£ • u>(C) 

13: for all a,b,c,d £ V such that Pl hcd is legal do 

14: for 2 <— 1 to A; do 

15: obtain G' from G by contracting the paths of P^ b d 

16: obtain w' from w by removing the ith objective 

17: P TSP MC-MaxATSP(G', w', k-l,e) 

18: for all H' G P TSP do 

19: form a Hamiltonian cycle from H' plus P^hcS a dd it to Ptsp 

20: form a Hamiltonian cycle from H' plus (u,v); add it to Ptsp 

Algorithm 3: Approximation algorithm for /c-criteria Max-ATSP. 



4.1 Multi-Criteria Max-ATSP 

We first present our algorithm for Max-ATSP (Algorithm 3) since it is a bit easier to analyze. 

First of all, we compute a (1 — e) approximate Pareto curve C of cycle covers. Then, for 
every cycle cover C £ C, we decide whether it is a light-weight cycle cover or a heavy-weight 
cycle cover (line 7). If C has only light-weight edges, we decompose it to obtain a collection 
P of paths. Then we add edges to P to obtain a Hamiltonian cycle H, which we then add to 
Ptsp- 

If C contains a heavy-weight edge, then there exists an edge e = (n, v) and an i with 
Wi(e) > afc • Wi(C). We pick one such edge. Then we iterate over all possible vertices a,b,c,d 
(including equalities and including u and v). We denote by P^bcd g ra ph with vertices u, 
v, a, b, c, d and edges (a,u), (u,b), (c,v), and (v,d). We call P^bcd ^ e 9 a ^ if it can be extended 
to a Hamiltonian cycle: P^bcd is legal if and only if it consists of one or two vertex-disjoint 
directed paths. Figure 2 shows the different possibilities. 

For every legal Pl bcd , we contract the paths as follows: We remove all outgoing edges of a 
and c, all incoming edges of b and d, and all edges incident to u or v. Then we identify a and 
b as well as c and d. If P^ b j consists of a single path, then we remove all vertices except the 
two endpoints of this path, and we identify these two endpoints. 

In this way, we obtain a slightly smaller instance G' . Then, for every i, we remove the ith 
objective to obtain w', and recurse on G' with only k — 1 objectives w' . This yields approximate 
Pareto curves P XSP of Hamiltonian cycles of G' . Now consider any H' € P-psp- We expand 
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(a) Two disjoint paths. (b) One path with an intermediate vertex. (c) One path including e. 

Figure 2: The three possibilities of -P^ bcrf - Symmetrically to (b), we also have a = d. Sym- 
metrically to (c), we also have v = a and u = d. 

the contracted paths to obtain H. Then we construct two tours: First, we just add Pl hcd to 
H, which yields a Hamiltonian cycle by construction. Second, we observe that no edge in H 
is incident to u or v. We add the edge (u, v) to H as well as some more edges such that we 
obtain a Hamiltonian cycle. We put the Hamiltonian cycles thus constructed into Ptsp- 

We have not yet discussed the success probability. Randomness is needed for computing 
the approximate Pareto curves of cycle covers and the recursive calls of MC-MaxATSP with 
k — 1 objectives. Let N be the size of the instance at hand, and let Pk(N, 1/e) is a polynomial 
that bounds the size of a (1 — e) approximate Pareto curve from above. Then we need at most 
iV 4 -pk(N, 1/e) recursive calls of MC-MaxATSP. In total, the number of calls of randomized 
algorithms is bounded by some polynomial <7fc(iV, 1/e). We amplify the success probabilities of 
these calls such that the probability is at least 1 — 2 - qk (N Tie) • Thus, the probability that one 
such call is not successful is at most qk(N, 1/e) ■ 2- qk (N,i/ £ ) ^ 1/2 by a union bound. Hence, 
the success probability of the algorithm is at least 1/2. 

Instead of Decompose, we can also use its randomized counterpart RandDecompose. We 
modify RandDecompose such that the running-time is guaranteed to be polynomial and 
that there is only a small probability that RandDecompose errs. Furthermore, we have to 
make the error probabilities of the cycle cover computation as well as the recursive calls of 
MC-MaxATSP slightly smaller to maintain an overall success probability of at least 1/2. 

Overall, the running-time of MC-MaxATSP is polynomial in the input size and 1/e, which 
can be seen by induction on k: We have a polynomial time approximation algorithm for k = 1. 
For k > 1, the approximate Pareto curve of cycle covers can be computed in polynomial 
time, yielding a polynomial number of cycle covers. All further computations can also be 
implemented to run in polynomial time since MC-MaxATSP for k — 1 runs in polynomial 
time by induction hypothesis. 

Theorem 4.1. MC-MaxATSP is a randomized jrjr[ — e approximation for multi- criteria 
Max-ATSP. Its running-time is polynomial in the input size and 1/e. 

Proof. We have already discussed the error probabilities and the running-time. Thus, it re- 
mains to consider the approximation ratio, and we can assume in the following, that all ran- 
domized computations are successful. We prove the theorem by induction on k. For k = 1, 
this follows since mono-criterion Max-ATSP can be approximated with a factor 2/3 > 1/2. 

Now assume that the theorem holds for k — 1. We have to prove that, for every Hamiltonian 
cycle H, there exists a Hamiltonian cycle H E Ptsp with w(H) > (jA^ — e) -w(H). Since every 
Hamiltonian cycle is in particular a cycle cover, there exists a C G C with w(C) > (l-e)-w(H). 
Now we distinguish two cases. 

The first case is that C consists solely of light-weight edges, i. e., w(e) < • w(C), then 
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Decompose returns a collection P of paths with w(P) > A^ -w(C) > (A^ -e) -w(H), which 

yields a Hamiltonian cycle H with w(H) > w(P) > {-Ai — e) ■ w(H) as claimed. 

The second case is that C contains at least one heavy-weight edge e = (u,v). Let (a, u), 
(u,b), (c,v), and (v,d) be the edges in H that are incident to u or v. (We may have some 
equalities among the vertices as shown in Figure 2.) Note that H does not necessarily contain 
the edge e. We consider the corresponding P^bcd an< ^ divide the second case into two subcases. 

The first subcase is that there exists a j G [k] with Wj{P^ b c d ) > -Ai ■ Wj(H), i. e., at least 
a jrjrf fraction of the jth objective is concentrated in P^ bcd - (We can have j = i, but this is 
not necessarily the case.) Let J C [k] be the set of such j. 

We fix one j & J arbitrarily and consider the graph G' obtained by removing the jth objective 
and contracting the paths (a,u,b) and (c, v,d). A fraction of 1 — A-^ = A-^ of the weight of 

H is left in G' with respect to all objectives but those in J. Thus, G' contains a Hamiltonian 
cycle H' with w e (H') > ^ • w e {H) for all I G [k] \ J. Since (k - l)-criteria Max-ATSP can 
be approximated with a factor of | — e by assumption, "P^sp contains a Hamiltonian cycle H' 
with we(H') > (| - e) ■ fejr ■ ^(#) > - e) ■ w t (H) for all i e [k] \ J. Together with 

Pabcd> w hich contributes enough weight to the objectives in J, we obtain a Hamiltonian cycle 
H with w(H) > (jAi — e) ■ w(H), which is as claimed. 

The second subcase is that Wj(Pabc d) — FFT ' w j(H) ^ OT ah j £ [k]. Thus, at least a fraction 
of -j^j of the weight of H is outside of P^ bcd - We consider the case with the ith objective 
removed. Then, with the same argument as in the first subcase, we obtain a Hamiltonian cycle 
H' of G' with we(H') > (-Ai — e) ■ Wf (H) for all I G [k] \ {i}. To obtain a Hamiltonian cycle of 
G, we take the edge e = (u, v) and connect its endpoints appropriately. (For instance, if a, b, c, d 
are distinct, then we add the path (a,u,v,d) and the edge (c,b).) This yields enough weight 
for the ith objective in order to obtain a Hamiltonian cycle H with w(H) > (^- — e) • w(H) 

since Wi(e) > ^ • w(C) > (^ - e) ■ w(H). □ 
4.2 Multi-Criteria Max-STSP 

MC-MaxATSP works of course also for undirected graphs, for which it achieves an approxi- 
mation ratio of — e. But we can do better for undirected graphs. 

Our algorithm MC-MaxSTSP for undirected graphs (Algorithm 4) starts by computing an 
approximate Pareto curve of cycle covers just as MC-MaxATSP did. Then we consider each 
cycle cover C separately. If C consists solely of light-weight edges, then we can decompose C 
using Decompose. If C contains one or more heavy-weight edges, then some more work has 
to be done than in the case of directed graphs. The reason is that we cannot simply contract 
paths - this would make the new graph G' (and the edge weights w') asymmetric. 

So assume that a cycle cover C G C contains a heavy-weight edge e = {u, v}. Let i G [k] be 
such that Wi(e) > Wi{C)/k. In a first attempt, we remove the ith objective to obtain w' . Then 
we set w'(f) = for all edges / incident to u or v. We recurse with k — 1 objectives on G with 
edge weights w'. This yields a tour H' on G. Now we remove all edges incident to u or v of H' 
and add new edges including e. In this way, we get enough weight with respect to objective i. 
Unfortunately, there is a problem if there is an objective j and an edge / incident to u or v such 
that / contains almost all weight with respect to Wj : We cannot guarantee that this edge / is 
included in H without further modifying H' . To cope with this problem, we do the following: 
In addition to u and v, we set the weight of all edges incident to the other vertex of / to 0. 
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Vtsp <- MC-MaxSTSP(G, w, k, e) 


input: undirected complete graph G = (V, E), k > 2, edge weights w : E — > N , e > 


output: approximate Pareto curve "Ptsp f° r fc-criteria Max-TSP 


1 


compute a (1 — e) approximate Pareto curve C of cycle covers 


2 


•p Tqp <_ 

' 1 or 


3 


if fc = 2 then 


4 


for all C G C do 


5 


P <— Decompose(C, fc) 


6 


add edges to P to form a Hamiltonian cycle .ff 


7 


add H to Ptsp 


8 


else 


9 


for all cycle covers C G C do 


10 


if iy(e) < w(C)/k for edges e G C then 


11 


P < — Decompose((7 w k) 


12 


add edges to P to form a Hamiltonian cycle 


13 


add H to Ptsp 


14 


else 


15 


let i € ffcl and e = {u,v\ £ C with wAe] > wAC)/k 


16 


for all £ G {0, . . . , 4k}, distinct xi, . . . ,x# G V \ {u, v}, and k G [A;] do 


17 


U <- {xi, . . . ,x^,u,w} 


18 


obtain w' from w by removing the jth objective 


19 


set w'( f) = for all edges / incident to U 


20 


pVgp MC-MaxSTSP(G, «/, fc - 1, e) 


21 


for all H G P^SP do 


22 


remove all edges / from H with / C U to obtain i?' 


23 


for all i7r/ such that i7' U Hjj is a Hamiltonian cycle do 


24 


add H' U F{/ to Ptsp 


Algorithm 4: Approximation algorithm for /c-criteria Max-STSP. 



Then we recurse. Unfortunately, there may be another objective j' that now causes problems. 
To solve the whole problem, we iterate over all £ = 0, . . . , 4/c and over all additional vertices 
xi, . . . , X£ 7^ u, v. Let U = {x\, . . . , X£, u, v}. We remove one objective i G [fe] to obtain w' , set 
the weight of all edges incident to U to 0, and recurse with k — 1 objectives. Although the 
time needed to do this is exponential in k, we maintain polynomial running-time for fixed k. 

As in the case of directed graphs, we can make the success probability of every randomized 
computation small enough to maintain a success probability of at least 1/2. 

The base case is now k = 2: In this case, every cycle cover possesses a 1/2 decomposition, 
and we do not have to care about heavy-weight edges. Overall, we obtain the following result. 

Theorem 4.2. MC-MaxSTSP is a randomized \ — e approximation for multi- criteria Max- 
STSP. Its running-time is polynomial in the input size and 1/e. 

Proof. We have already dealt with error probabilities and running-time. Thus, we can assume 
that all randomized computations are successful in the following. What remains to be analyzed 
is the approximation ratio. As in the proof of Theorem 4.1, the proof is by induction on k. 

The base case is k = 2. Let H be an arbitrary Hamiltonian cycle. Then there is a C G C 
with w(C) > (l-e)-w(H). From C, we obtain a Hamiltonian cycle H with w(H) > ^-w(C) > 
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(\ — e) • w{H) by decomposition and Lemma 3.3. 

Let us analyze MC-MaxSTSP for k > 3 objectives. By the induction hypothesis, we can 
assume that MC-MaxSTSP is a ^rry — £ approximation for (k — l)-criteria Max-STSP. Let H 
be any Hamiltonian cycle. We have to show that "Ptsp contains a Hamiltonian cycle H with 
w(H)>(l-e)-w(H). 

There is a C G "Ptsp with w{C) > (1 — e) ■ w(H). We have to distinguish two cases. First, 
if C consists solely of light-weight edges, i. e., w{e) < w{C)/k for all e, then we obtain a 
Hamiltonian cycle H from C with w{H) > w{C)/k > (| — e) ■ w(H). 

Second, let e G C and i G [k] with Wi(e) > Wi(C)/k. We construct sets I C [k], X C H, and 
(/ C 1/ in phases as follows (we do not actually construct these sets, but only need them for 
the analysis): Initially, I = X = and U = {u, v}. In every phase, we consider the set X' of 
all edges of H that have exactly one endpoint in U. We always have < 4 by construction. 
Let I' = {j G [k] | j I,Wj(X') > Wj(H)/k}. If I' is empty, then we are done. Otherwise, 
add I' to /, add X' to X, and add all new endpoints of vertices in X' to U. We add at least 
one element to I in every phase. Thus, \X\ < 4k and \U\ < 4k + 2 since |/| < k. 

Let w m = w(X), and let w 9 = Y2feH-\fnu\=i w ^ ^e ^ e we ight, of edges of H that have 
exactly one endpoint in U. Let w ont = w(H) — w m — w 9 . By construction, we have w 9 < 1/k 
for all j ^ I. Otherwise, we would have added more edges to X. 

We distinguish two subcases. The first subcase is that 7 = 0. Then w m = and w 9 < 1/k. 
Consider the set "P^gp an< ^ the edge weights w' used to obtain it. We have w'j(H) = w° nt > 
(^r-) • Wj(H) for j i. By the induction hypothesis, there is an H G T^-jsp w ith 

* (^i -') ■ (M) ' »<*> £ (l " E ) • »<*> 

for j 7^ i. We remove all edges incident to u or v to obtain i?'. Since the weight of all these 
edges has been set to 0, we have w'(H') = w'(H). There exists a set Hq such that e G H$ and 
H' U is a Hamiltonian cycle. For this cycle, which is in Vtsp, we have 

Wi (H' U H 9 ) > Wi(e) > Wi{C)/k > Q - ^ ■ «;(£) 

and, for j 7^ i, 

wj(fr' U H ) > w'j(H) > Q - ■ «;(£). 

The second subcase is that I is not empty. Let j G I, and let {/. We consider V^ p . Let 
t« m , and u> out be as defined above. By the induction hypothesis, the set P T g p contains 
a Hamiltonian cycle cover H with w' e (H) > (j^ - e) ■ w^ ut for £ + j. We remove all edges 
incident to U from H to obtain H' with w'(H') = w'{H). By construction H'UX is a collection 
of paths. We add edges to X to obtain Hjj such that H' U Hjj is a Hamiltonian cycle. Let us 
estimate the weight of H' U R\j. For all i G I, we have we(H' U ilf/) > we(Hu) > we(H)/k. 
For all £ ^ J, we have 

u Ft,) > + < > - • (<"?"' + <4") > fjrry - E ) ■ - '"I) 

which completes the proof. □ 
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5 Remarks 



The analysis of the approximation ratios of our algorithms is essentially optimal: Our approach 
can at best lead to approximation ratios of for some c E Z. The reason is as follows: 
Assume that (k — l)-criteria Max-TSP can be approximated with a factor of Tfc. If we have a 
/c-criteria instance, we have to set the threshold for heavy-weight edges somewhere. Assume 
for the moment that this threshold be arbitrary. Then the ratio for fe-criteria Max-TSP is 
minjafc, (1 — • Tfc_i}. Choosing = T ^ k ~\i maximizes this ratio. Thus, if Tfc_i = l/T for 

some T, then < T Tfc ~^. 1 = tTj- We conclude that the denominator of the approximation 
ratio increases by at least 1 if we go from k — 1 to k. 

For undirected graphs, we have obtained a ratio of roughly 1/k, which is optimal since 
a% = 1/2 implies c > 0. Similarly, for directed graphs, we have a ratio of ^rpy, which is also 
optimal since a 1 ^ = 1/S implies c > 1. 

Due to the existence of 0(1/ log /^-decompositions, we conjecture that both fc-criteria Max- 
STSP and /c-criteria Max-ATSP can in fact be approximated with factors of 0(1/ log k). This, 
however, requires a different approach or at least a new technique for heavy-weight edges. 
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